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1. Introduction 

In recent years, it has become an established fact that neutrinos, though relatively light, 
are massive. Since the first experimental evidence of neutrino oscillations until today an 
enormous progress has been made in determining the low-energy properties of neutrinos, 
such as mass squared differences and mixing. The existence of neutrino masses poses some 
fundamental theoretical challenges, such as understanding why the neutrino mass is so 
much smaller than the masses of the other fermions. An elegant and attractive solution to 
this problem is given by the seesaw mechanism [1—9], which explains the smallness of the 
neutrino mass through the existence of very heavy particles (usually right-handed Majorana 
neutrinos or Higgs triplets) , the mass scale of which could be related to that of Grand Uni- 
fication. In addition, the seesaw mechanism provides a natural explanation of the baryon 
asymmetry of the Universe through the baryogenesis via leptogenesis mechanism [10] (for 
recent reviews, see refs. [11-13]). However, the large mass scale of the seesaw particles 
jeopardizes the hopes of testing this mechanism in the laboratory and hence reduces its 
predictivity. 

In the present work, we consider the seesaw mechanism in a class of left-right symmetric 
models in which the intermediate states with both right-handed neutrinos (type I) and 
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heavy triplet scalars (type II) contributions to the light neutrino mass matrix m u are 
naturally present. We focus on a special case with a discrete left-right symmetry, in which 
type I and type II seesaw contributions contain the same triplet Yukawa coupling /. This 
case has much fewer parameters than the most general one and is therefore more predictive. 
After integrating out the heavy particles, the light neutrino mass matrix is given by 



v 2 

m u = fv L yf~ 1 y T , (l.l) 



VR 



where / is the triplet Majorana-type Yukawa coupling, y is the Dirac-type Yukawa coupling 
of neutrinos and v, Vl, and vr are vacuum expectation values (VEVs). The first term in 
eq. fll.lD is the type II contribution, while the second term is the type I contribution from 
the original seesaw scenario. In the case when y is a complex symmetric matrix, it was 
shown in ref. [14] that if the light neutrino mass matrix m v , the VEVs, and the Dirac-type 
Yukawa coupling matrix y are known, the seesaw formula (1.1) can be inverted analytically 
to find the triplet Yukawa coupling matrix /. Since the seesaw equation is non- linear in 
/, one can expect multiple solutions, and indeed an eight-fold of allowed solutions is found 
[14]. As the mass matrix of heavy right-handed Majorana neutrinos is given by Mr = fvR, 
this also implies an eight-fold ambiguity for this mass matrix. For given Dirac-type Yukawa 
coupling matrix y and VEVs, all eight solutions for / result in exactly the same mass matrix 
of light neutrinos m„, and thus, the seesaw relation by itself does not allow one to select 
the true solution among the possible ones. One therefore has to invoke some additional 
information and/or selection criteria. The present work is an attempt in this direction. 

One possibility to discriminate among the eight allowed solutions for / is to introduce 
a notion of naturalness. For example, for certain ranges of the VEVs and certain solutions, 
a very special triplet coupling matrix / might be needed, in the sense that marginally 
different / would lead to significantly different low-energy phenomenology. We consider 
such a situation unnatural; the degree of tuning that is required in the right-handed sector 
to obtain the observed neutrino phenomenology will be quantified and the corresponding 
selection criterion for / discussed in section ||. 

Another possibility to discriminate among the allowed solutions is to constrain them 
by the phenomenology of the right-handed neutrinos. Since the right-handed sector of the 
theory is not directly accessible to laboratory experiments, cosmological benchmarks turn 
out to be the most promising tool. Namely, we will classify the solutions according to 
their ability to lead to successful baryogenesis via leptogenesis. This will be discussed in 
section f|, before we draw our conclusions in section || 

Recently, leptogenesis in a class of models with the left-right symmetric seesaw mech- 
anism has been considered in a similar framework in ref. [15]. We compare our results with 
those in ref. [15] in sections || and [|. 
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2. The model and the inversion formula 

In this section, we introduce our framework and set up the notation. In the basis where the 
mass matrix of charged leptons is diagonal, the light neutrino mass matrix can be written 
as 

m u = (P t U PM ns Pu)* mf* g (Pi £/pmns Pu) ] , (2.1) 

where m^ iag = diag(mi, 7712, W3) is the diagonal matrix of neutrino masses, t/pMNS is 
the leptonic mixing matrix which depends on three mixing angles and a Dirac-type CP- 
violating phase, and Pi and P u are diagonal matrices of phase factors, which in general 
contain five independent complex phases. 

The neutrino masses mi, m,2, and 711,3 can be expressed through the lightest neutrino 
mass mo and the two mass squared differences Am^ and Am|j . In our numerical calcula- 
tions, we will use the current best-fit values of the parameters defining the neutrino mass 
matrix [16-18]: 

Am^ ~ 7.9 x 1(T 5 eV 2 , Am^ ~ ±2.6 x 1(T 3 eV 2 , (2.2) 

012 ~ 33.2° , # 23 ~ 45° . (2.3) 

For the mixing angle #13, only the upper limit #13 < 11.5° exists. Unless explicitly stated 
otherwise, we will use the value #13 = in our analysis. 

We will be assuming that the Dirac-type Yukawa coupling matrix of neutrinos y coin- 
cides with that of the up- type quarks y u . This is a natural choice in the light of quark- lepton 
symmetry and grand unified theories (GUTs) [19-21]. Actually, this relation is unlikely to 
hold exactly, since, in the GUT framework, it would also imply that the Yukawa couplings 
of the down- type quarks and charged leptons coincide, yd = yu in contradiction with ex- 
periment. GUT models that modify this relation usually also modify the relation between 
the up-type and neutrino Yukawa matrices [22, 23]. However, most of the qualitative re- 
sults in the present work depend only on the fact that the eigenvalues of y are hierarchical. 
Whenever a result relies on the assumption y = y u , we will comment explicitly on this 
issue. Following ref. [14], we will also assume y to be symmetric. In this case, the two 
VEVs (vl and vr), the sign of Am^, and the mass scale of the light neutrinos are the only 
free parameters (ignoring for the moment the CP-violating phases, which will be discussed 
in section |3|) . 

Our choice of the Dirac-type Yukawa coupling matrix implies that it can be written as 

y = Pd uEkm p u yt iag Pu UcKM Pd , (2.4) 
where the eigenvalues of yt ia,g are 

y£ iag = diag(4.2 x 10~ 6 , 1.75 x 1(T 3 , 0.7) , (2.5) 
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and we use the standard parameterization of the CKM matrix Uqkm [24] with 

q 12 ~ 13.0°, 9 q 13 ~ 0.2°, 0f 3 ~ 2.2°, 5 9 ~ 1.05 . (2.6) 
The values in eqs. (p.5|) and (|2.6|) are evaluated at the GUT scale, following ref. [15]. The 



matrices P u and Pd in eq. (2.4) are diagonal matrices of phase factors. The phases in the 
four matrices Pi, P u , P u , and Pd are partially redundant. For example, by a redefinition 
of the fields, the three phases of Pi can be moved into Pd, so that we are left with the two 
usual Major ana phases and the Dirac phase in the low-energy sector, while five additional 
Majorana phases and one Dirac phase reside in y and can only affect high-energy processes 
such as leptogenesis. Even though these phases can marginally influence the stability of 
the seesaw solutions, we set the high-energy phases to zero in the first part of our work 
and consider them only in the part where leptogenesis is discussed. 

In order to invert the seesaw formula, it is useful to introduce the following dimensionful 
quantities: 

9 = vif, M = —^o ' ( 2 - 7 ) 

vl f 



with the VEV v ~ 174 GeV, so that eq. (|l.l| ) turns into 



m v = g--yg l y T . (2.8) 
A* 

This convention has the advantage that the matrix g will only depend on \x and not on the 
two VEVs, vl and vr, separately. It will turn out that the baryon asymmetry produced via 
leptogenesis depends only on this combination of VEVs, so that, besides the CP-violating 
phases, we are left with two parameters only, the quantity [i and the lightest neutrino mass 
tuq . The hierarchy of the light neutrino masses can be considered as an additional discrete 
parameter. 

In the following, we give a short description of the seesaw inversion formula from 
refs. [14, 25] in the case of three lepton generations and when y is a complex symmetric 
matrix. In the basis where y is diagonal, the seesaw equation for g reduces to the following 
system of six coupled non- linear equations for its matrix elements g^ : 

fj,G\gij - (m u )ij] = yiyjGij . (2.9) 

Here we use the notation 

1 3 

G = detg, Gy = - ^ikl^jmn9km9ln ■ (2-10) 

k,l,m,n=l 

It was found in ref. [14] that in the case when y is symmetric, for every solution g there 
exists another solution g which is related to g by the duality transformation g = m u — g. 
For g, eq. ( ^S| ) reads 

fj.G[gij - (m u )ij] = -fiGgij = yiyjGij (2.11) 
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with G = det g. The system of equations in eq. ( pT9|) can now be solved by making use of the 
following procedure. First, we introduce the rescaled matrices g' = g/X 1 ^ 3 , m' v = m^/A 1 / 3 , 
and y' = y/X 1 ^ 3 , where A is to be determined from the equation G"(A) = detg'(X) = 1. 
Then, using the equation for the dual quantities g' , one can linearize the system of equations 
for g'-. Next, this system can be solved and one obtains the following solution for g: 



9ij 



where 



A 2 [(A 2 - Y 2 ) 2 - Y 2 Xdetm u + Y i S]{m v ) lj + A(A 4 - Y )Aij - Y 2 X 2 (X 2 + Y 2 )S, 



i.i 



(A 2 - Y 2 ) 3 - Y 2 X 2 (X 2 - Y 2 )S - 2Y 2 X 3 det m„ 



Y A 



(2/12/22/3) 



/r 



k,l=l 



{m v ) 2 



Ikl 



VkVl 



A 



ViVjMij 



ft 



k,l=l 



kl 



VkVl 



(2.12) 
(2.13) 

(2.14) 



with Mij 



2 e iki£jm,n{'rni>)krn{iTi'v)in- In terms of the original (non-rescaled) quantities, one 
has G(A) = det g(X) = A, which yields an eighth order equation for A. Using the duality 
property, one can reduce it to a pair of fourth order equations. Substituting the solutions 
for A into eq. ( 2.12| ) gives eight solutions for g$j. In general, for n lepton generations the 
number of solutions is 2 n [14]. 

The matrix structure of the solutions of the seesaw equation was studied in some 
detail in ref. [25]. In the present work, we will rather focus on the eigenvalues of the 
matrices g, the corresponding mixing parameters, stability properties of the solutions, and 
the implications for leptogenesis. 



3. Stability analysis 

Since the neutrino Dirac-type Yukawa coupling matrix in our framework is given by the 
up-type quark mass matrix, the inversion formula of the previous section can be used to 
determine the eight allowed structures of the triplet coupling matrix / = g/v^ for given low- 
energy neutrino mass matrix m u and the parameters vl, vr, and itlq. Our stability analysis 
is based on the assumption that the Dirac-type coupling matrix y and the Majorana-type 
coupling matrix / are a priori independent (for a discussion of the situations when this 
is not the case, see section 5 of ref. [25]). We pose the question of whether the resulting 
low-energy phenomenology is stable under small changes in /. Since the inversion formula 
in general yields eight valid solutions, the mass matrix m v and the corresponding Majorana 
coupling matrix / are in a l-to-8 correspondence. It is still a reasonable question to ask if 
for the measured m v some of the predicted / have to be very special, so that a fine-tuning 
is required and a small modification of their elements may lead to a large change in (m u )ij. 
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The measure we use to quantify the stability property of the solutions is the following: 

det/ 



Q 



l/n 



(3.1) 



det m u 

The real coefficients f k and m; determine the matrices / and m„ according to 

/ = ^(/fc + */fc+iv)Tfc, (3.2) 

(3.3) 



in,, 



^2(m k +im k+N )T k , 
k 



where T k , k G [l,iV] with iV = n(n + l)/2, form a basis of complex symmetric n x n 
matrices. For this basis, we choose the normalization 



tr {TjT k ) = 5 t 



ik 



(3.4) 



The resulting stability measure Q does not depend on the chosen basis. This can be easily 
seen in the following way. Consider another basis T' k satisfying eq. ( |3.4| ). The two bases 
are then connected via a unitary transformation TL = U k i T\. The coefficients in the old 
and new bases are determined as 



fk 



Re 



tr (Tl f) 



fkA 



^=Re[tr(^/)l, f kA 



■N 



■N 



Im 
Im 



tr(r' f /) 



and hence, are related by an orthogonal transformation 



f'a = Y,O ab f b , a,be[l,2N], 0=(^ 



ReU ImU 
-ImU ReU, 



(3.5) 
(3.6) 

(3.7) 



which leaves the measure in eq. ( |3.1| ) invariant 1 . 

Many interesting properties of the seesaw inversion formula appear already in the 
one-flavor case. The solutions g are then given by 



m,, 



ml y2 
4 + n 



(3.8) 



and our stability measure simplifies to 



Q = f 37 log \m v \ = g ^- log \m v 
df dg 



1 + 



(3.9) 



In the following, we will discuss the qualitative behavior of the solutions / in various regions 
of the parameter space and its implications for the stability of these solutions. 



Note that the stability issue was also discussed in ref. [15] where a different stability criterion, con- 
straining only the element /33, was introduced. 
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Figure 1: An example of our labeling convention for the solution '— 



10 



10' 



> 10 

a io 

,40 



10° 

io 6 



10 



4x10" 



2x10" 



io 14 io 16 io 18 io 20 io 22 io 24 io 26 





III 












/ / 




/ / 




1 / 




/ 




/ 

*s 



v /v 

R' L 



io 14 io 16 io 18 io 20 io 22 io 24 io 26 

v /v 

R L 



Figure 2: The right-handed neutrino masses mjq i and mixing parameters as functions of vr/vl 
for the solution 1 '. Normal mass hierarchy, mo = 0.001 eV. 



3.1 Large /i regime 

In the regime of large /x, 



/i > — 2 



the two solutions in the one-flavor case are given by 

,2 



9 



and (7 



(3.10) 



(3.11) 



In this regime, the solutions are purely type I or type II dominated. In the three-flavor 
case, the eight solutions follow from the eight corresponding choices for the eigenvalues and 
we will label these solutions according to their limiting behavior at large \i as '— ' or '+' in 
the case of type I or type II dominance (starting with the largest eigenvalue in the small fj, 
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Figure 3: Same as in fig. H. but for the solution 
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Figure 4: The right-handed neutrino masses mjVi and mixing parameters as functions of vr/vl 
for the solution 'H h'. Inverted mass hierarchy, m = 0.001 eV. 
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Figure 5: Same as in fig. 0, but for the solution 



regime). This notation agrees with the one used in ref. [15]. Our convention is illustrated 
in fig. [l] using the solution ' — h +' as an example. 

From eq. ( |3.9[ ) one can observe that in the large fx regime of the one-flavor case, both 
solutions for g are characterized by the stability measure Q ~ 1, which is a very stable 
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situation. Note that for the three-flavor case, no fine-tuning corresponds to Q ~ 10. 
However, for several flavors and hierarchical y, there is in general an instability related to 
mixing that will be discussed in the next subsection. 

3.2 Hierarchy induced large mixing 

For simplicity, we start with a discussion of the two-flavor case in the pure type I seesaw 
framework. By hierarchy induced large mixing we mean the following: Suppose that y has 
a hierarchical structure 

(3.12) 

while, in contrast to this, the low-energy neutrino mass has a rather mild or even no hier- 
archy. Then, the corresponding matrix g is necessarily characterized by the hierarchy that 
is the squared hierarchy of y. Indeed, introducing a unitary matrix U (6) that diagonalizes 
g, one finds 

g = --ymZ 1 y = U\e)gU*(d) (3.13) 
I* 

with the diagonal matrix 

»~(o?)- <3 ' 14) 

and, in addition, mixing has to be small, i.e. 9 ~ e. This was already observed in refs. [26- 
28] and suggested as a possible mechanism for generating large mixing angles in the light 
neutrino mass matrix out of small mixing angles in the right-handed and Dirac sectors. 
However, in our context, this is not a desirable situation, since it would require a fine- 
tuning between the Dirac and Majorana Yukawa couplings, i.e. between the sectors that 
we have assumed to be unrelated. In terms of stability, this would lead to large values of 
Q. In addition, the large hierarchy among the elements of the Dirac-type Yukawa coupling 
matrix y would induce a huge hierarchy among the elements of g, leading in general to an 
extremely small mixing in the right-handed neutrino sector, which may preclude successful 
leptogenesis. 

The above consideration was based on the type I seesaw formula, and hence, is not 
fully applicable to our framework. Still, it applies to the solutions dominated by type I 
seesaw. Figure ^ shows the one out of the eight solutions that is fully dominated by the 

type I term in the large fi regime and is labeled as ' '. As a measure of mixing, we 

consider the parameters U{ which are related to the off-diagonal elements of the unitary 
matrix U diagonalizing g as follows 2 : 

"l = ^(l^l2| 2 + |C/2l| 2 ), ^2 = ^(l^l3| 2 + |C/3l| 2 ), n! = i(|£/ 23 | 2 + |[/ 32 | 2 ). (3.15) 



2 Recall that we work in the basis where the matrix y is diagonal. 
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These parameters, along with the masses of right-handed neutrinos, are plotted for several 
solutions in figs. |2||| 

For the solution ' ', mixing is small in the large /U regime, as can be seen from 

fig. For the other seven solutions, this does not hold in general, as can be seen e.g. in 
fig. ||. However, even in the general case, one feature seems to be universal: If the matrix 
elements of g exhibit a strong hierarchy, then the mixing in the right-handed sector is 
suppressed, which leads to the necessity of fine-tuning between the Dirac and Majorana 
sectors and related instabilities. This also explains why the two solutions '+ H — ' and 

'H ' are very unstable with almost equal stability measure Q. The strong hierarchy 

between the largest and smallest right-handed masses leads to large instabilities, while the 
behavior of the third mass is rather irrelevant. 

3.3 Small /x regime 

When /i is small in the sense that 

V , , 

H < -V , (3.16) 
in the one-flavor case, one finds the following limiting behavior for g: 

9-±-^ + ^ + 0(v^), (3.17) 



1,1 ' oo (3.18) 



For the stability measure, eq. Q3.9| ) gives 

g dm u 



Q 



m v dg 



in this limit, and therefore a very unstable situation. This had to be expected, since there 
is an almost exact cancellation between the type I and type II contributions to m v in the 
seesaw formula in this regime. In the multi-flavor case, there is an additional instability in 
the small fi limit which stems from the fact that mixing in g is suppressed by the hierarchy 
in y. This can be illustrated by the two-flavor case, in which the four solutions are of the 
form 

a = -^y 1/2 Py 1/2 (3.19) 

with P of the form 

(cos ex sin (y. \ 
+ C7(v7t). (3.20) 
sin a — cos a J 

For the first pair of solutions, mixing vanishes in the limit fj, — * 0, while for the second 
pair, mixing in g is suppressed by the hierarchy in y. A similar argument applies to the 
three-flavor case and can be observed in our numerical results. For example, this behavior 
can be seen in figs. |2| and ||] which display two out of the eight solutions for the normal 
mass hierarchy and mo = 0.001 eV. 
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Figure 6: The stability measure Q as a function of vr/vl for mo = 0.001 eV. The left (right) 
panel corresponds to the normal (inverted) neutrino mass hierarchy. 



3.4 Numerical results 



Figures || and show the stability measure Q for small and large and normal/inverted 
mass hierarchy. For small mo, the transition from the large fi to the small fi regime appears 
for larger values of /i, in accordance with eqs. ( |3.10| ) and ( 3.16| ). In all four scenarios, the 
solutions are unstable in the regime of small fi, which is due to the cancellation between 
type I and type II contributions to the mass matrix of light neutrinos. In addition, the 
solutions where the smallest eigenvalue is dominated by type I seesaw in the large /i regime 
('db =b — '), are unstable for large [i as well, since the lightest right-handed mass stays below 
10 6 GeV in this limit and this generally leads to a large spread in the eigenvalues and to 
instabilities, as explained in the previous sections. Examples of the eigenvalues in these 
cases are given in fig. Analogously, the stability measure of the solutions '± — h' increases 
for vr/vl > 10 20 , since the smallest right-handed neutrino mass approaches its asymptotic 
value of about 10 9 GeV, as can be seen in figs. || and || A similar effect appears for the 
solution ' — h +' at values vr/vl > 10 24 . The purely type II dominated solution ('+ + +') 
is the most stable one in almost all the cases. If one allows for a tuning at a percent level, 
Q < 10 3 , then the stability analysis favors the two solutions 'db + +' with vr/vl > 10 18 
and the two solutions '± — h' with vr/vl — 10 20 . 

It should be noted that the qualitative behavior of the stability measure Q depends 
mostly on the eigenvalues of the Yukawa coupling matrix y and the neutrino mass scale 
mo- On the other hand, the mixing in y and additional CP-violating Majorana phases 
influence the results only marginally. 
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Figure 7: Same as in fig. g, but for m Q — 0.1 eV. 

4. Leptogenesis 

In this section, we present our analysis of leptogenesis and its implications for the discrim- 
ination among the eight allowed solutions for g. Our analysis is based on the results of 
refs. [29, 30]. 

Assuming that the lightest of the right-handed neutrinos is separated from the other 
two as well as from the Higgs triplets by a large mass gap, the baryon asymmetry arising 
from leptogenesis can be written as 

T] B = — =r]e Nl . (4.1) 

n 7 

The observed value of the baryon asymmetry is t]b = (6.1 ± 0.2) x 10~ 10 [31]. In eq. (f4.1|), 
r\ is the so-called efficiency factor that takes into account the initial density of right-handed 
neutrinos, the deviation from equilibrium in their decay and washout effects, while 
denotes the lepton asymmetry produced in the decay of the lightest right-handed neutrino. 
For the decay of the ith right-handed neutrino, it is defined as 

e _ T(N i ^lH)-T(N i ^lH*) 

Ni T(Ni^lH) + T(Ni^lH*)' { ' } 



If the two lightest right-handed neutrinos have similar masses, eq. (4J.) is generalized 

to 

VB = m eJVi + m £7V 2 • ( 4 -3) 
The coefficients r\i mostly depend on the effective neutrino masses, defined as 

v 2 {$y)u .x 

to,- = . 4.4 

2m Ni K ' 
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Figure 8: The contributions 771 and 772 to the baryon-to-photon ratio from the decays of the two 
lightest right-handed neutrinos versus the ratio of their masses m^/m^. Left panel: rh\ = 7712 = 
10~ 3 eV, right panel: rh\ = 10~ 3 eV, m 2 = 10~ 2 eV. 



Here and below, the hat indicates that the matrices are evaluated in the basis where the 
triplet Yukawa coupling matrix g is diagonal with real and positive eigenvalues. In the case 
of quasi-degenerate right-handed neutrinos, — mjv 3 , and nearly coinciding effective 
masses rh\ and r/7.2, an order-of-magnitude estimate of the washout coefficients gives [32] 

1 /10 " 3eV ^ (4.5) 



' 200 V rhi 

However, deviations from the condition rh\ ~ rfi2 can lead to large corrections to this 
estimate. In particular, a large effective mass 7712 reduces the coefficient 771 close to the mass 
degeneracy point, as is shown in fig. || The results in ref. [32] have been obtained for rather 
light and quasi-degenerate right-handed neutrinos, ~ ~ 1 TeV. For hierarchical 
right-handed neutrino masses <C ttin 2 and the mass scale under consideration in the 
present case, mjv a ~ 10 8 GeV, one finds 

/ in-3 e y\ 

7/1 = 1.45 x 10~ 2 — z , r, 2 ~ , (4.6) 



mi 

and we will employ these values in the following. This result and fig. ^ have been obtained 
by solving the Boltzmann equations as suggested in ref. [32] and assuming thermal initial 
abundance of right-handed neutrinos. 

With the washout factors rji at hand, the determination of the baryon asymmetry 
requires only the knowledge of the CP-violating decay asymmetries of the right-handed 
neutrinos e^- In the case when the low-energy limit of the theory is the Standard Model, 
tN 1 is given by [30] 

= e J Nl + e$ , (4.7) 
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II 3 , MB'SlHu] 



1 — z In 



* + l 



(4.9) 



and analogous formulas hold for e^ 2 . Here z = m\/m 2 Ni , and Xj is denned as the ratio of 
the squared right-handed neutrino masses: 



9JI 
9 2 n 



(4.10) 



In the following, we discuss only the limit of a very heavy SU(2)l Higgs triplet, z 
so that 

u . 3 A Jm[(yt^*) u ] 



OO, 



(4.11) 



167r Jiir " (yty)ii 

In the limit of a strong hierarchy in the right-handed sector, xj 3> 1, the first contribution 
in eq. (4/7) can be rewritten as 



1 ^lm[(y*y)y\ 



8tt 



ii 



3 „ lm[($ yg 
5ii — 



■W)u] 



16vr 



ii 



so that 



+ e 



n 

Ni 



3 , Im[(yt m^) n ] 

I6vr (y f y)n 



(4.12) 



(4.13) 



However, even in this limit, this approximation can lead to large deviations from the exact 
result of eqs. (f4T^)-( |4~9| ). Consider e.g. the regime of small /x, where type I and type II 
seesaw contributions almost cancel each other in the expression for the light neutrino mass 
matrix. In this case, even a small correction to the coefficient of the asymmetry ejy leads 
to an incomplete cancellation and to large errors in the approximation of eq. ( [4.13 ). This 
effect is also partially present at intermediate values of [i. In addition, close to the mass 
degeneracy (xj ~ 1), a resonant feature is expected in ejy , which can lead to successful 
leptogenesis even at a TeV scale [32]. This is demonstrated in fig. pi, where the asymmetries 
6 jVi an d €n 2 produced in the decays of the two lightest right-handed neutrinos and the 
corresponding effective mass parameters rh\ and rfi2 are plotted. The results show sizable 



deviations from the approximation (4.13), even outside the resonant enhancement region. 



The corresponding baryon-to-photon ratio is shown in fig. 1C . In addition, this figure shows 
the baryon-to-photon ratio in the case of non- vanishing #13 and the Dirac-type leptonic CP- 
violating phase (5cp = 30°. The resonant behavior is less distinct for larger values of #13, 
which can be traced back to the fact that the two lightest right-handed neutrinos never 
become exactly degenerate in mass in this case. On the other hand, the Dirac-type phase 
constitutes an additional source of CP violation in the case of non- vanishing #13, leading 
to an enhancement of ejvi below the mass degeneracy point for smaller values of #13, and 
thus, widening the vr/vl region where successful leptogenesis is possible (see the dashed 



curve in fig. |10| ). 
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Figure 9: The upper (lower) panels show the effective neutrino mass rhi/eY (rfi2/eV) and the 

asymmetry ejvi ( £ n 2 ) as functions of vr/vl for the solution 'H \-\ The dashed curves in the right 

panels correspond to the approximation in eq. (4.13), while the solid curves represent the exact 
result. The step-like behavior of fh\ and TO2 is due to the level crossing. Inverted mass hierarchy, 
toq = 0.001 eV. 



Thus, we find that viable leptogenesis is possible in this scenario if the ratio of the 
VEVs is close to vr/vl ~ (1 4- 2) x 10 19 . Note that leptogenesis in the case of the left-right 
symmetric seesaw mechanism was previously considered in a similar framework in ref. [15]. 
For the specific choice of the parameters made there, the washout processes were found to 
be too strong to allow successful leptogenesis. However, for our choice of the parameters 
with the inverted mass hierarchy in the light neutrino sector, the drop in the effective mass 
mi below the level crossing point of the two lightest right-handed neutrinos resolves this 
issue. We notice that the use of the exact formulas (4.7-4.9) rather than the approximation 
(4.13) is essential in this region. 

It should be also noted that a similar effect of incomplete cancellation can appear if 
the mass of the Higgs triplet is of the same order as the mass of the lightest right-handed 
neutrino. In this case, the asymmetry e J is modified and the cancellation between type 
I and type II contributions is incomplete as well, which in the small and intermediate 
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Figure 10: The baryon-to-photon ratio t/b from the decay of the lightest right-handed neutrino 

for the solution H h'. The same parameters as in fig. |[ except that the dashed and dotted curves 

correspond to nonzero #13 and Sep = 30°. The shaded area corresponds to values of t]b below the 
observed value. 



\jl regimes can enhance the produced lepton asymmetry by several orders of magnitude 
compared to the approximation in eq. ( 4.13j ). 



With the parameters of fig. |10|, the lightest right-handed neutrino has a mass of order 
m N x — 5 x 10 9 GeV, as can be seen in fig. ||. Since thermal leptogenesis requires a reheating 
temperature T > Mjsr 1: this can potentially lead to a tension with bounds coming from 
gravitino cosmology in supersymmetric theories, namely T < (10 7 -j- 10 10 ) GeV [33]. Thus, 
this possibility imposes constraints which are similar to those in the usual pure type I 
seesaw scenario. 

Another difference from the standard leptogenesis scenario is the appearance of the 
phases contained in P u , Pi, P u , and Pd in the neutrino mass matrix m u and in the Dirac 
Yukawa coupling matrix y, which up to now have been set to zero in our discussion. Due 
to these phases and an interplay between type I and type II contributions to the neutrino 
mass matrix, leptogenesis is possible, in principle, even in the case of one leptonic flavor, 
as will be demonstrated below. This case is quite similar to the framework with three left- 
handed neutrinos and one right-handed neutrino discussed in ref. [34] (see also ref. [25]). 
In the following, we will present some analytic results for the left-right symmetric one- and 
two-flavor cases, before presenting numerical results for the three-flavor case. 

In the one-flavor case, the light neutrino mass is given by 

V 2 

m v = g - — (4.14) 
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and the lepton asymmetry produced in the decay of the heavy right-handed neutrino is 

J_lm[r 2 m,] 

32vr m V ; 

Once again, the hat indicates that y and m„ are in the basis where g is real and positive. 
It turns out that the most interesting regime is given by large values of fx and a relative 
phase of 7r/4 between m u and y. In this case, only the solution dominated by the type II 
term is relevant, since the type I contribution to y* 2 m u is real and cannot generate any CP 
asymmetry. Thus, we obtain 

g ~ m„, iriN = m^fiv 2 , (4.16) 

and 

ft _tefV = J 2 !L i (4 . 17) 

IrriN lm v \x 
3 2 3 m u m N 

6 = 16^ = 1^;^' 2 (418) 

77b = 1.7 x 10- 6 eV ^ = 1.7 x lO" 6 eV ^| . (4.19) 

M \y\ v 

Thus, it is possible to reproduce the observed baryon asymmetry e.g. with the values 

|y| = 10- 4 , m„ = 0.1eV, n = 6.0 x 10~ 5 eV" 2 , (4.20) 

which leads to 

m = 8.3 x 10~ 4 eV, m N = 1.8 x 10 8 GeV . (4.21) 

The situation, however, is more complicated in scenarios with more than one lepton 
flavor. For instance, mixing could give large contributions to mi, thereby enhancing the 
washout. On the other hand, it can also lead to additional sources of CP violation, which 
might improve the prospects for successful leptogenesis in realistic models with several 
flavors. Consider, for example, the situation when the third right-handed neutrino is much 
heavier than the other two and the mixing with the third flavor in the right-handed sector 
is suppressed. A novel aspect of this effective two-flavor case is that large mixing and 
resonant amplification of the lepton asymmetries due to the level crossing of right-handed 
neutrinos can enhance leptogenesis. These effects are similar to those discussed above in 
the full three-flavor framework. We will study the regime with a large hierarchy between 
the two lightest right-handed neutrinos, which allows a simple analytic approach. As a 
toy example, we consider the following scenario: We assume maximal mixing in the light 
neutrino sector and one complex phase in Pi, which can be moved into the Yukawa coupling 
matrix y by rephasing the electron neutrino field. Thus, the neutrino mass matrix is taken 
to have the form 

(e 2m m e tK Sm\ 
■ e (4-22) 
e lK dm m 
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with 5m <C m. The parameters m and 5m can be determined from the mass of the lightest 
active neutrino mo and Am^: 

Ami . . 

m ~ mo, om ~ — — — . (4.23) 
4m 

Numerical analysis indicates that the most interesting region in the parameter space cor- 
responds to the situation when the smaller eigenvalue of g is in the large \x regime, while 
the larger eigenvalue is in the small \x regime, i.e. 

4 4«,«'4, (4.24) 
m z m z 

and we will assume this to hold in the present example. In this case, two solutions for g 
are, to first order in A, given by the ansatz 3 

g = U ] \ , _\U\U=\ .. 4.25 
±4 + t \ -\e*+ 1 



with 



VP 2 



X = zp Y-t- , (4.26) 

V2 

sin(0 + k) ~ =F sin(2K)^— , (4.27) 

m^/JI 



and thus, we find 



rh, = yl + vlX2 = yJ±^ , (4.28) 
2m/i 2m[i 
3 

e^, = \sm(2(fr + 2k) m 5m 2 /j, + sm(4K)myi] . (4.29) 

327rmi L J 

The second term in ejvi essentially coincides with the corresponding expression in the one- 
flavor case. Hence, in this case, it is possible to generate a sufficient lepton asymmetry in 
exactly the same way as in the one-flavor case as long as the contribution from mixing to 
mi does not lead to a strong washout. The latter condition reads 

5m 2 ^ (Arn^) 2 ^ in _ 3 



2m 2,2ml 



<10" d eV, (4.30) 



which is easily satisfied if mo > 1CP 3 eV. It is interesting to note that for k = 7r/8 and 
quasi-degenerate neutrino masses, the obtained asymmetry ejv a saturates the upper limit 
obtained in ref. [30]. 

But even in the case k ~ 7r/4, when the second term in the expression for in 



eq. ( 4.29D is suppressed, the first term can lead to viable leptogenesis. The corresponding 



3 The other two solutions do not lead to successful leptogenesis. 
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contribution to t\b takes its largest value when Sm 2 = y\/ Li, so that eqs. ( [4.28 ) and ( 4.2S| ) 
become 

mi = (4.31) 
mil 

3 _ 

= T^-y^Vi 1 - ( 4 - 32 ) 

lbvT 

In this case, tjb is smaller than it is in the one- flavor case only by a factor 

yi i 




2^ <"- 33 > 
It should be noted that the baryon asymmetry increases with the parameter ii, so that, 



depending on the Yukawa couplings, saturation of the upper limit on /i, in eq. (424) might 
be necessary, which can lead to deviations from our analytic results. 

Thus, in the two-flavor case, two different sources of leptogenesis exist: The first source 
is similar to that in the one-flavor case, which is related to the type II seesaw term and 
is sensitive to the high-energy CP-violating phases, while the second source results from 
mixing effects and has no analogue in the one-generation case. 

In the three-flavor framework, sources of both types are, in general, present as well, 



but mixing with the third flavor can further increase rh\. Figure |ll| shows the baryon-to- 
photon ratio r/s when an additional phase is attributed to the electron neutrino, as in the 
two-flavor example of eq. (|4.22| ). We choose the phase k = 7r/4 (k = tt/8), so that the 



source similar to the first (second) term in eq. ( 4.29 ) gives the largest contribution to the 
baryon asymmetry. Our numerical results indicate that, similarly to the two-flavor case, the 
upper bound on the decay asymmetry found in ref. [30] can be saturated. The mass of the 
lightest right-handed neutrino that is required to reproduce the observed baryon asymmetry 
is niN 1 > 1-4 x 10 9 GeV (m^ > 2.5 x 10 8 GeV). These bounds can be relaxed by choosing 
Yukawa couplings different from those of the up-type quarks. With an appropriate choice, 
the results for the four solutions of the type '± ± +' agree with the analytic predictions 
of the two-flavor analysis presented in this section. Notice that the results in the two- 
flavor case in eqs. ( [4.28 ) and ( 4.29| ) do not depend on y2 as long as the constraint ( 4.24| j 



is fulfilled. Likewise, we observe in the numerical analysis of the three-flavor case that in 
this limit leptogenesis is not very sensitive to the two largest eigenvalues y<i and y%. This 
is, however, a consequence of the fact that the mixing in the 1-3 sector of the Dirac-type 
Yukawa coupling y is small in our framework according to eq. fl2.6|) . If this mixing is sizable, 
$13 ^ 5°, and depending on the other parameters determining the Yukawa coupling y and 
the neutrino mixing matrix C/pMNS) leptogenesis might be suppressed, mainly due to a large 
contribution to the effective mass parameter rh\ from the eigenvalue 7/3 and the resulting 
increased washout. 
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Figure 11: The baryon-to-photon ratio r\B with an additional complex phase 7r/8 or 7r/4 attributed 

to the electron neutrino for the solution ' Y\ The shaded area corresponds to values of t\b below 

the observed value. Inverted mass hierarchy, too = 0.1 eV. 



Thus, we conclude that successful leptogenesis is possible for four out of the eight 
solutions provided that the value of the electron-type Majorana phase is in an appropriate 
range. For the other four solutions, leptogenesis is not viable, as was first pointed out 
in ref. [15]. The reason for this is that, as long as the Dirac Yukawa coupling matrix is 
chosen to coincide with that of the up-type quarks, the mass of the lightest right-handed 
neutrino never exceeds 10 6 GeV and no level crossings occur. We note that in the left-right 
symmetric case with type I+II seesaw mechanism the bounds on the mass of the lightest 
right-handed neutrino can be slightly relaxed compared to those in the pure type I case 
which, for right-handed neutrinos with thermal initial abundance and hierarchical masses, 
requires m Nl > 5 x 10 8 GeV [35-37]. 



5. Summary and conclusions 





± + + 


±- + 


±±- 


Stability 


v R /v L > 10 18 


vr/v l ^ 10 20 


disfavored 


Leptogenesis 


vr/vl > 10 18 


vr/v l > 10 18 


excluded 


Gravitinos 


vr/v l < 10 21 


unconstrained 


unconstrained 



Table 1: The allowed regions of the parameter vr/vl for the eight different types of solutions. 

We have analyzed the left-right symmetric type I+II seesaw mechanism with a hierar- 
chical Dirac mass term motivated by GUTs. It was previously shown that a reconstruction 
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of the mass matrix of heavy right-handed neutrinos in this framework produces eight so- 
lutions which result in exactly the same low-energy phenomenology. Our goal was to 
discriminate among these solutions using their stability properties and leptogenesis as ad- 
ditional criteria. As a measure of the stability, we have chosen the parameter Q which 
quantifies the degree of fine-tuning necessary to obtain a given mass matrix of light neutri- 
nos and was defined in eq. ( |3.1| ). For three lepton generations, no fine-tuning corresponds 
to Q ~ 10. We have selected the value Q = 10 3 , which corresponds to a fine-tuning at 
the percent level, as a maximal allowed value. The leptogenesis criterion we used was the 
ability of a given solution to reproduce the observed baryon asymmetry of the Universe. 

Our results complement the results of the leptogenesis analysis performed in ref. [15] 
in the following aspects. In the case without additional Majorana phases, we obtain, in 
accordance with ref. [15], that a sizable decay asymmetry is possible close to the mass 
degeneracy of the two lightest right-handed neutrinos. However, while for the specific 
parameters used in ref. [15] the washout is too large to allow viable leptogenesis, we find 
that assuming the inverted mass hierarchy for the light neutrinos resolves the problem, as 
shown in fig. ^. Similarly, in the cases with additional CP-violating Majorana phases we 
found that for certain solutions the choice of the parameters made in ref. [15] leads either 
to a strong washout (solutions '± — (-')) or to a violation of the gravitino bound (solutions 
'± + +'). In section we presented a systematic study showing that those problems can 
be solved for the four solutions '± ± +' if the value of the of electron- type Majorana phase 
is in the appropriate range. In particular, the upper bound on the decay asymmetry for 
the type I+II seesaw model found in ref. [30] can be saturated for a certain choice of the 
parameters. This is illustrated by the analytic results for the two-flavor case in eqs. ( 4. 2g ) 
and ( 4.29| ) and the numerical results for the three- flavor case in fig. 11. We would like to 
emphasize that if the Dirac-type Yukawa coupling matrix y is characterized by hierarchical 
eigenvalues and rather small mixing, successful leptogenesis is quite a generic feature of 
the left-right symmetric seesaw models. 

Our findings are summarized in tab. |l[ One can observe that the stability criterion 
disfavors the four solutions of the type '± ± — ' and restricts the solutions of the type 
'± — h' to the region of the parameter space where vr/vl — 10 20 . The remaining two 
solutions of the type '± + +' are stable, provided that vr/vl > 10 18 . We found that 
successful leptogenesis is possible for the four solution of the type '± ± +' as long as 
v r/ v l ^ 10 18 . This possibility requires the existence of additional Majorana- type phases 
which are absent in the pure type I seesaw framework. Further constraints come from 
the potentially dangerous overproduction of gravitinos in supersymmetric theories, giving 
rise to an upper bound on the lightest right-handed neutrino mass. For our choice of the 
Yukawa couplings, y = y u , only the solutions of the type '± + +' are affected by this 
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constraint, which leads to the requirement vr/vl < 10 21 . For the other six solutions, the 
smallest right-handed neutrino mass is always below 10 10 GeV, so that these solutions 
are not constrained by this criterion. In the cases when the middle eigenvalue of y is 
chosen to be significantly larger than the one in our framework, yi > 10" 2 , the constraint 
vr/vl < 10 21 would also apply to the two solutions of the form '± — \-\ On the other 
hand, a very small middle eigenvalue, y2 < 5 x 10~ 4 , would render leptogenesis impossible 
for these two solutions, since the decay asymmetry would be too small due to the small 
mass of the lightest right-handed neutrino. 

Thus, we have shown, within the chosen framework, that the stability and leptogenesis 
criteria partially lift the eight-fold degeneracy among the solutions for the mass matrix of 
heavy right-handed neutrinos in the left-right symmetric type I+II seesaw. 
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